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ABSTRACT 

We perform a series of direct V-body calculations to investigate the effect of resid¬ 
ual gas expulsion from the gas-embedded progenitors of present-day globular clusters 
(GCs) on the stellar mass function (MF). Our models start either tidally filling or un¬ 
derfilling, and either with or without primordial mass segregation. We cover 100 Myr 
of the evolution of modeled clusters and show that the expulsion of residual gas from 
initially mass-segregated clusters leads to a significantly shallower slope of the stel¬ 
lar MF in the low- (to ^ O. 5 OM 0 ) and intermediate-mass (~ 0.50 — 0.85Mq) regime. 
Therefore, the imprint of residual gas expulsion and primordial mass segregation might 
be visible in the present-day MF. We find that the strength of the external tidal field, 
as an essential parameter, influences the degree of flattening, such that a primordially 
mass-segregated tidally-filling cluster with ry/rt ^ 0.1 shows a strongly depleted MF 
in the intermediate stellar mass range. Therefore, the shape of the present-day stellar 
MF in this mass range probes the birth place of clusters in the Galactic environment. 
We furthermore find that this flattening agrees with the observed correlation between 
the concentration of a cluster and its MF slope, as found by de March! et ah. We show 
that if the expansion through the residual gas expulsion in primordial mass segregated 
clusters is the reason for this correlation then GGs most probably formed in strongly 
fluctuating local tidal fields in the early proto-Milky Way potential, supporting the 
recent conclusion by Marks & Kroupa. 

Key words: globular clusters- initial mass segregation - methods: N-body simula¬ 
tions 


1 INTRODUCTION 

The study of dense stellar systems, which mainly refers to 
young massive star clusters, globular clusters and ultra¬ 
compact dwarf galaxies, has played an important role in 
developing our knowledge about the universe. Since most 
stars in the Galactic disc may originate in embedded star 
clusters, these systems can be viewed as the fundamental 
building blocks of galaxies to understand the origins of the 
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properties of the galactic stellar population, such as the 
galaxy-wide s tellar mass function dWeidner &: Kroupall200^ : 
iKroupalboOSl l. 

The internal properties of star clusters can undergo sig¬ 
nificant changes at birth and also during the course of cluster 
evolution. Hence, it is important to specify to what amount 
the current properties of star clusters are imprinted at the 
beginning of the evolution and early formation processes and 
to what amount they are the result of long-term evolutionary 
processes (Vesperini 2010). In fact, star formation provides 
the initial conditions for the long-term evolution of stellar 
clusters, while the follow up dynamics of stars in clusters 
can inform on the conditions that have dominated in the 
environment in which the clusters have formed. 
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The stars of star clusters are born more or less simulta¬ 
neously within the dense cores of the same progenitor giant 
molecular cloud Q. Stars within an embedded cluster form 
following an initial mass function (IMF). The shape of the 
stellar IMF is one of the main uncertainties in our under¬ 
standing of star formation, i.e., the conversion of gas into 
stars over cosmic time. It is usually assumed that the IMF in 
external galaxies is the same as in the disk of the Milky Way, 
where most studies of resolved stellar populations showed 
the IMF to be invariant (Bastian et al. 2010, Kroupa et al. 
2013). The IMF is also observationally found to be invariant 
for nearly all dense systems in which its shape could be in¬ 
ferred and which formed with a star formation density below 
a threshold near Q.IMq/{ yrpc?) and is hence referred to as 
the "standard” or "canonical” IMF. This poses a problem 
for star formation theories which predict a dependence on 
the environment where star formation takes place (Kroupa 
et al 2013). 

In the course of time, the MF of stars in clusters evolves 
at early stages and also during the cluster secular evolu¬ 
tion as the stellar members are removed from the system 
through stellar and dynamical evolution (Vesperini & Heg- 
gie 1997, Baumgardt & Makino 2003). Simulations show 
that mass segregation and the preferential loss of low-mass 
stars which leads to a shallower mass function should only 
happen after a cluster has gone into two-body relaxation- 
driven core-collapse which is connected to the shrinkage of 
the core size rc (Baumgardt & Makino 2003). So one would 
expect a correlation between MF-slope, a {dN/dm oc m~°‘, 
a = 2.3 being the Salpeter index), and concentration param¬ 
eter (c = log(rt/rc)), in the sense that clusters with large 
values of c are depleted in low-mass stars (fig. 4 in Leigh et 
al. 2013). Here, dN is the number of stars in the mass in¬ 
terval m to m -I- dm and rt, Vc are the tidal and core radius, 
respectively. 

But, according to obser vational data by D e Marchi, 
Paresce, & Pulone (2007), and iPaust et al.l (l2010l l for Galac¬ 
tic globular clusters the value of a, in the mass range 0.3-0.8 
Mq , correlates surprisingly with their concentration param¬ 
eter, such that the less concentrated clusters are typically 
but not exclusively more strongly depleted in low mass stars. 
This trend can not be understood to be the result of long¬ 
term evolution alone and at least some initial correlation 
between the initial concentration and the total cluster mass 
is required (Leigh et al 2013). 

The ’De Marchi relation’ could be interpreted as a first 
evidence for a non-universality of the IMF at low masses, 
but Baumgardt et al (2008a) showed that such a correlation 
can be understood if clusters are born initially mass segre¬ 
gated with a canonical IMF, but are filling their tidal radii. 
The difficulty of this scenario is that the cluster would have 
to form with knowledge of its mass-dependent tidal radius at 
different locations in a galaxy. It is unclear how star forma¬ 
tion can arrange this. Marks et al. (2008), however, demon¬ 
strated, using N-body integrations of young, gas embedded 
compact clusters, that the loss of the primordial residual- 
gas from initially mass-segregated clusters in a range of tidal 


^ However, many observations suggest that massive clusters ap¬ 
pear to be composed of (at least) two distinct stellar populations 
(e.g., Piotto 2010) 


field strenghts starting with a canonical IMF and containing 
unresolved binaries reasonably reproduces the relation be¬ 
tween the slope of the actual mass function and the concen¬ 
tration of globular clusters within the observational limits. 

The flattening of the MF-slope driven by a violent early 
phase of gas-expulsion of an embedded cluster with primor¬ 
dial mass segregation can also be a possible explanation for 
the observed flattened mass function of some remote halo 
globular clusters like for example, Palomar 4 and Palomar 
14 (Zonoozi et al. 2011, 2014) 

In this paper, our aim is to consider the effect of primor¬ 
dial residual-gas expulsion on the early evolution of globular 
clusters (covering the first 100 Myr of their evolution) and 
it’s influence on the stellar MF-slope and degree of mass seg¬ 
regation for a wide range of stellar masses, varying the star 
formation efficiency, and the strength of the external tidal 
field. Contrary to the otherwise identical models studied by 
Marks et al. (2008), we here use a fully realistic distribution 
of stellar masses. Our paper is organized as follows: In Sec. 2 
we discuss the star cluster reaction to residual-gas expulsion. 
In Sec. 3 we describe the initial set-up of the N-body models 
and in Sec. 4 we present the main results of the simulations. 
Finally Sec. 5 consists of the conclusions and a discussion. 


2 THE REACTION OF STAR CLUSTERS TO 

RESIDUAL-GAS EXPULSION 

The fraction of gas that is converted into stars within the 
volume of the embedded cluster is defined as the star for¬ 
mation efficiency (SFE), 

^ _ Mecl 

^ecl T Algas 

where M^ci is the stellar mass of the embedded cluster and 
Mgas is the residual gas within the cluster. A large SFE leads 
to more stars and therefore also to more massive stars with 
stronger winds and radiation and there is less gas to remove. 
The observationally determined SFE is usually smaller than 
e ~ 40% (Lada & Lada 2003) so most of the mass of the 
progenitor cloud remains in the form of residual gas to be 
expelled. 

In young, still gas-embedded star clusters, the left-over 
gas from star formation is lost over a time which is deter¬ 
mined by feedback from massive stars through UV radiation 
and massive stellar winds from OB stars or supernova explo¬ 
sions. Thus, star clusters will become super-virial and their 
further dynamical evolution will be strongly affected by the 
gas loss. 

Recent observations indicate that proto-stars form in 
molecular clouds as complex and filamentary substructures. 
However, as shown in Banerjee & Kroupa (2015), if a young 
star cluster has to form out of them, they must fall in the po¬ 
tential well of the molecular cloud and form a near-spherical 
and nearly virialized cluster within a few dynamical times, 
which is short (or similar) compared to the typical gas dis¬ 
persal time that we take here (« 0.6 Myr). Observational ev¬ 
idence for such a monolithic formation of individual massive 
cluster rather than a hierarchical formation of successively 
merging of small-clusters is provided by the > 5 x 10® M©, 
and < 24 pc sized molecular cloud formed in the Anten¬ 
nae merging galaxies by Johnson et al. (2015). Hence, our 
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initial conditions, that comprise virialized embedded clus¬ 
ters, are plausible. As suggested by some recent theoretical 
stndies involving hydrodynamic calculations, e.g.. Dale et 
al. (2015, 2012), the locations at which the protostars (or 
sink particles) are formed, the star formation efficiency is 
nearly 100% and so the removal of gas is less efficient in 
further dynamics, arguing against gas-expulsion as being a 
relevant physical process that governs the long-term evoln- 
tion of a cluster. This is likely to be a result of including 
only radiative feedback and ignoring magnetic fields as well 
as performing simulations in which the sink particles are 
significantly more massive than individual stars. However, 
according to the high resolution hydrodynamic calculations 
by Machida & Maatsnmoto (2012) and Bate et al. (2014), 
and direct observation (Lada & Lada 2003) only small frac¬ 
tion of gas convert into stars. 

Based on the hyd r odyna mical simulations without feed¬ 
back, iKruiissen et al.l (l2012tl also claimed that the SFE is 
higher in the regions where stars have formed, and hence 
the long-term evolution of star clusters is only weakly af¬ 
fected by gas removal. They argued that this is because 
of the low gas fractions within the pre-cluster cloud core, 
caused by the accretion of gas on to the stars. In this case 
it however needs to be explained how the outer parts of the 
gas cloud are expelled without heating and pushing out the 
inner gas contents of the pre-cluster cloud core. It also is 
difficult to explain ultracompact HII regions if gas can sim¬ 
ply be accreted onto the stars despite their luminosity. In 
the short-lived UCHII phase, the gas is heated significantly 
(up to K or more) in the immediate surrounding of the 
stars (within 0.05 pc), which then leads to an outflow of the 
residual gas on larger scales. 

However, _ in the simulations by 

IPale. Ercolano. fc Bonndl (l2012l l including the influ¬ 
ence of photoionizing radiation from 0-type stars it is 
shown that the global SFE is between 20-30% and agrees 
with what is observed in star-forming regions and embedded 
clusters which suggests that the star formation efficiency 
can at most be e « 30% (see, e.g., Lada & Lada 2003), 
although on the scale of whole molecular clouds the star 
formation efficiency is at most a few percent (Murray 
2011). High-resolution mapping in star-forming regions has 
uncovered that outflows, even from low-mass proto-stars, 
are very powerful in destroying and clearing cloud material 
away from the star-forming region on sub-ps scales (e.g., 
Bally et al. 2011, 2012, McGroarty et al. 2004, Frank et al. 
2014). 

In addition, the dynamical state of the stars at the 
onset of gas expulsion might b e a key factor in determin¬ 
ing the effects of gas expulsion llVerschueren fc Davidl[l989l : 
lGoodwinll2009l L For example, by means of N-body simula¬ 
tions of highly non-equilibrium star clusters in fixed back¬ 
ground potential (as a simple model of gas potential). Smith 
et al. (2011) found that the long-term dynamical evolution 
of star clusters is affected by the initial stellar distribution 
and velocity dispersion far more than by gas expulsion. They 
have shown that initially cool clusters are more able to sur¬ 
vive gas expulsion, weakening the dependency of survival to 
gas expulsion on the initial SFE compared to a virialised 
initial dynamical state of the clusters. Such an ansatz how¬ 
ever requires all stars to have formed in some volume syn- 
chronically on a time-scale much shorter than the dynamical 


time-scale of this volume, which appears to be an artihcial 
initial state. Further discussion of this issue is available in 
Kroupa (2005). Observations of UCHII regions, of outflows 
and magnetohydrodynamical simulations of star formation 
show that the SFE does not surpass about 20-30% on a sub- 
pc scale because matter is blown out in outflows launched 
by the magnetic fields. This questions the above mentioned 
results according to which gas-removal is not important. 

The main argument here is that actual proto-star scale 
high-resolution MHD calculations (e.g., Machida & Maat- 
sumoto 2012, Bate et al. 2014) show that the interplay be¬ 
tween gas accretion onto the protostar and magnetically 
driven outflow results in star formation efficiencies up to 
« 30% in the neighborhood of proto-stars, which is a well 
sub-pc scale. This would imply the overall SFE in a dense gas 
clump (sub-pc to pc scale) would be less than 30%, includ¬ 
ing its innermost regions. Such proto-star scale calculations 
are the only available ones that can track the mass accretion 
and outflow (driven by radiation and magnetic field) from in¬ 
dividual proto-stars. Other more massive SPH simulations, 
that intend to form star clusters, either exclude all feedback 
or include only partial feedback. The SFE from these calcu¬ 
lations can never be reliable. The Dale et al. (2015) calcu¬ 
lations that go up to globular cluster masses had their sink 
particle representing a sub-cluster itself, i.e., their sink par¬ 
ticles themselves are of pc scale. Of course, with a 30% SFE 
all over the clump (the SFE can, in practice, be smaller in 
the outer regions), gas expulsion should have a significant 
impact on the cluster. This is the initial star-to-gas ratio 
that we plausibly adopt in our computations. 

We assume that the gas expulsion starts after the em¬ 
bedded phase at a certain time, td, which is set equal to 
0.6 Myr and corresponds to the UCHII phase. After this 
time, it is assumed that the mass of gas decreases expo¬ 
nentially on the gas-expulsion characteristic time scale tge, 
over which the natal gas is removed from the cluster. So the 
time-varying total mass in gas is given by 

M,as{t) = Mgas{Q) (2) 

The impact of the early gas expulsion on the early dy¬ 
namical evolution of star clusters has been investigated by 
Baumgardt & Kroupa (2007) using a large set of N-body in¬ 
tegration. They showed that the survival rate and final prop¬ 
erties of star clusters are strongly influenced by the details 
of the expulsion of the residual-gas that did not form stars, 
depending on the parameter values for the gas-expulsion 
time scale, tge, the star formation efficiency, e, and the 
ratio of the initial half-mass radius to tidal radius. They 
also found that if the tidal field is weak {rn/rt ^ 0.05, ru 
being the half-mass radius) and the gas is removed slowly 
{tge tcross, tcross being the crossing time defined here 
as tcross = 2rhlo, where a is the 3D velocity dispersion in 
the cluster), star clusters can survive even if the SFE is as 
low as 10%. It is shown here that the gas-free stellar system 
(after gas expulsion) expands violently and loses a fraction 
of its stars depending on its initial mass and concentration 
if e « 1/3 and tge is the thermal gas removal timescale. 
We take tge = ’ where 10 km/s is the approximate 

sound velocity of the ionized gas. Studying the dynamical ef¬ 
fect of gas dispersal by adopting a time-varying background 
potential mimicking the influence of the gas with the above 
values for td tge and e is used by Banerjee & Kroupa (2013, 
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2014, 2015) to explain well observed very young massive 
clusters such as R136 and the NGC 3603, and the ONC is 
likewise well accounted for using these values as a currently 
expanding population leading to a Pleiades-like open cluster 
in ~ lOOMyr (Kroupa et al. 2001). The conformity of these 
values to account for observed clusters of masses 10®, 10"*, 
and 10® M© suggests a common physical process which acts 
in their emergence from their embedded stage. 


3 DESCRIPTION OF THE MODELS 

The effective dynamical influence of gas-removal from 
embedded clusters is included as modeling the gas 
as an additive time-varying potential and accord- 
ingly modifying the equat i ons o f motion of the stars 
dKroupa. Aarseth. fc Hurl^ I 2001 I : iBaumgardt fc Kroupal 
l200'it) . I.e.. gas was not simulated directly, instead we 
mimic the essential dynamical effects of the gas-expulsion 
process by applying a diluting, spherically-symmetric exter¬ 
nal gravitational potential to a model cluster. Specifically, 
we use the potential of the spherically-symmetric, time- 
varying mass distribution as described by Eq. 2. Comparing 
computations treating the gas with the SPH method, Geyer 
& Burkert (2001) have demonstrated that such an analytical 
approach is physically realistic. 

The models presented here assume the pre-gas expul¬ 
sion stellar plus gas configuration of the cluster models to 
be in virial equilibrium. This is a simplifying assumption 
which is necessary to make the models computable, because 
treating the gas component with a magneto-hydrodynamical 
(MHD) and radiation-transfer code including stellar feed¬ 
back (i.e., radiation and matter out flows) to the star form¬ 
ing gas is computationally impossible even for a small cluster 
of mass larger than about 100 Mq (Klessen et al. 1998; Bate 
& Bonnell 2004; Bate 2009; Girichidis et al. 2011, 2012; Bate 
2012). In reality the earliest phases of embedded cluster as¬ 
sembly are given by the simultaneous inward flow of gas, 
driven by the deepening potential well of the embryonic em¬ 
bedded cluster, and by the nearly-simultaneous outflow of 
gas driven by the increasing stellar feedback and magnetic 
fields generated in the proto-stellar accretion disks. 

We use the state-of-the-art collisional N-body inte¬ 
gration code NBODY6 (Aarseth 2003) on GPU comput¬ 
ers to perform direct N-body simulations. NBODY6 uses 
a 4th-order Hermite integration scheme and an individual 
timestep algorithm to follow the orbits of cluster members 
and invokes regularization schemes to deal with the inter¬ 
nal evolution of small-N subsystems. NBODY6 includes the 
tidal effects of the Galaxy, and stellar and binary evolu¬ 
tion using the dynamical integrated SSE/BSE routines de¬ 
veloped by Hurley et al. (2000, 2002). It includes the time- 
varying gas potential to model the gas-expulsion process. 
The initial condition of all clusters are set up using the pub¬ 
licly available code McLuSTEF0(Kupper et al. 2011). 

In this work, all models are located at a Galactocentric 
distance of Rg = 8.5 kpc orbiting on a circular orbit with 
circular velocity of Vg = 220 km/s. The effect of the exter¬ 
nal tidal field is adopted to be the same as the Milky Way 

https://github.com/ahwkuepper/mcluster.git 


potential including bulge, disk and a Galactic NEW halo po¬ 
tential (Navarro et al. 1996, 1997) which can be viewed as 
a rough proxy for the log-normal isothermal phantom dark 
matter halo expected in Milgromian dynamics (Famaey & 
McGaugh 2012; Kroupa 2015). We calculate several mod¬ 
els with different initial cluster half-mass radius {vh =1, 5, 
7 and 10 pc) to cover two different regimes: tidally-filling 
{rhjrt ^ 0.1) and tidally-underfilling (rhjrt ^ 0.1) clusters. 
The Th «1 pc clusters are consistent with the radius - mass 
relatio n for embedded clusters inferred bv iMarks fc Krounal 
(l2012tl . We emphasize that the adoption of the galactic 
model here merely serves the purpose of defining a com¬ 
putable model. In reality, the Galaxy would not have ex¬ 
isted when the GGs formed and the proto-Galactic tidal field 
would have been evolving locally rapidly (Marks & Kroupa 
2010). 

The number of stars and initial cluster mass in all com¬ 
puted models are set to N=170000 and M = 10® Mq, re¬ 
spectively. The initial mass function of stars is adopted to 
be a canonical Kroupa IMF which consists of two power laws 
with slope q= 1.3 for stars with masses, m, between 0.08 and 
0.5Mq and a slope a=2.3 for more massive stars. Stellar evo¬ 
lution is modelled according to the routines by Hurley et al. 
(2000). The range of stellar masses was chosen to be from 
0.08 to IOOMq and the metallicity is Z= 0.001 dex (which is 
a typical value for the young halo GGs). The evolution time 
for all models is 100 Myr. Frank et al. (2012) estimate an 
age of 11 ± 1 Gyr for Pal 4 (as an example of the typical ob¬ 
ject we study) by adopting [Fe/H]= —1.41 TO. 17dex for the 
metal licity from the best-fitting isochrone of iDotter et al.l 

(l2008ll . 

For initially segregated models, the mass segregation 
was set up by using the segregati o n rou tine described in 
IBaumgardt. De Marchi. fc Kroup^ ll2008lf . The degree of 
segregation (parameterized by S) can be chosen to be be¬ 
tween 0 and 1, where S = 0 means no segregation and S = 1 
refers to full segregation, i.e. the most massive star sitting in 
the lowest energy orbit and the second-most massive stars 
sitting in the second-lowest energy orbit, and so on. This 
routine allows to maintain the desired density profile when 
increasing the degree of mass segregation while also making 
sure that the cluster is in virial equilibrium. The primor¬ 
dial segregation for all models are set to S=0.9, except one 
model which is set to S=0 for comparison. We do not have 
primordial binaries in our simulated clusters, although bi¬ 
naries created via three-body interactions are automatically 
included. Including a significant fraction of primordial bina¬ 
ries is computationally impossible for such massive clusters 
as moddeled here (see Leigh et al 2015 for details). 

Concerning the assumption of primordial mass segre¬ 
gation, recently, Parker et al. (2015) used the final state of 
SPH simulations of star formation which include feedback 
from photoionization and stellar winds as the initial condi¬ 
tion of their N-body simulations and found that after 10 Myr 
of evolution, primordial mass segregation is less common in 
the simulations when star formation occurs under the influ¬ 
ence of feedback. They found that the ten most massive stars 
are mass segregated, rather than full mass segregation of the 
whole cluster. These simulations concern low-mass clusters 
with very simplified treatment of feedback, such that it is 
unclear how indicative they are. For example, the outflows 
launched even by low-mass stars appear to be destructive on 
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Table 1. The initial and the final parameters of the N-body runs after t=100 Myr evolution. Column 1 gives 
the model name, in which the first two digit numbers after ‘rh’ denote the initial half-mass radius in parsec, 
the last two digit numbers denote the percentage of the star formation efficiency (sfe), and the degree of mass 
segregation is written after the letter ’S’. The initial mass of all models is M = IO^Mq. Column 2 gives the 
initial 3D half-mass radius. Columns 3, 4 and 5 contain the adopted star formation efficiency, gas expulsion 
timescale and initial degree of segregation, respectively. Column 6 is the initial filling factor, where rto = 54.3 pc 
for all models. Note that the gas mass is not included in the calculation of the tidal radius. The final slope of the 
mass function after 100 Myr of evolution in the mass range 0.50 Si ^ 0.85 inside the tidal radius is presented 
in column 7, the initial canonical value being 2.3. The following columns give the final mass, the 3D half-mass 
radius of the simulated star clusters after 100 Myr of evolution, the final concentration c = log^Q (rt/rc), as well 
as the two-body relaxation time of the star clusters after 100 Myr of evolution. 


Model 

rho 

[pc] 

£ 

tge 

[Myr] 

s 

{rh/n)o 

a2 

[0.50 - 0.85]Mo 

[Mq] 

rhf 

[pc] 

c 

Trh 

[Gyr] 

(1) 

(2) 

( 3 ) 

( 4 ) 

( 5 ) 

(6) 

( 7 ) 

(8) 

( 9 ) 

(10) 

(11) 

rhl0-S0-sfe25 

10 

0.25 

1 

0.0 

0.18 

2.16 

36741 

22.3 

1.38 

13.14 

rhl-S0.9-sfe33 

1 

0.33 

0.1 

0.9 

0.02 

2.27 

67392 

4.4 

2.40 

1.81 

rhl-S0.9-sfe25 

1 

0.25 

0.1 

0.9 

0.02 

2.21 

65016 

5.2 

2.10 

1.50 

rh3-S0.9-sfe25 

3 

0.25 

0.1 

0.9 

0.05 

1.85 

34930 

12.9 

1.55 

4.40 

rh5-S0.9-sfe25 

5 

0.25 

1 

0.9 

0.09 

1.96 

54274 

18.6 

1.61 

9.36 

rh7-S0.9-sfe25 

7 

0.25 

1 

0.9 

0.13 

1.45 

41716 

21.6 

1.34 

11.02 

rhl0-S0.9-sfe33 

10 

0.33 

1 

0.9 

0.18 

1.08 

31486 

22.6 

1.195 

11.81 

rhl0-S0.9-sfe25 

10 

0.25 

1 

0.9 

0.18 

0.81 

20060 

21.5 

1.178 

10.08 

rhl0-S0.9-NoGE 

10 

1.0 

- 

0.9 

0.18 

2.28 

66038 

16.2 

1.70 

8.84 

rhl-SO-NoGE 

1 

1.0 

- 

0.0 

0.02 

2.29 

73827 

1.9 

2.65 

10.10 

rh3-S0-NoGE 

3 

1.0 

- 

0.0 

0.05 

2.27 

74379 

4.3 

2.44 

12.40 

rh5-S0-NoGE 

5 

1.0 

- 

0.0 

0.09 

2.30 

72668 

7.0 

2.42 

8.01 

rhlO-SO-NoGE 

10 

1.0 

- 

0.0 

0.18 

2.22 

71039 

12.7 

1.84 

7.50 

rh20-S0-NoGE 

20 

1.0 

- 

0.0 

0.36 

2.26 

55029 

22.8 

1.67 

6.30 

rh30-S0-NoGE 

30 

1.0 

- 

0.0 

0.54 

2.34 

10864 

35.0 

0.86 

3.40 


a cloud scale. Observations of some young clusters such as 
ONC (Allison et al 2009), NGC3603 (Pang et al 2013) and 
Trumpler 14 (Sana 2010) have shown these clusters to be 
mass segregated down to around the 4*^, 20*^, and 6*^ most 
massive stars, respectively, rather than primordial mass seg¬ 
regation of the whole cluster. Thus, while the complete mass 
segregation in young clusters has not been observed yet, in 
this paper, we aim to address possible initial conditions of 
star-burst clusters in a totally different mass regime a Hub¬ 
ble time ago in low-metallicity gas. 

In our calculations, we assume that the SFE does not 
depend on the position inside the cluster, so both stellar and 
gaseous components follow the same density distribution ini¬ 
tially, which is given by a Plummer model with similar ra¬ 
dius. Using a Plummer model is a reasonable physical choice 
because it is the simplest solution of the collisionless Boltz- 
man equation and also adequately describes the simplest 
stellar-dynamical objects, globular clusters, known to exist 
(Plummer 1911) We assume different SFEs: e : 25%, 33%. 
The gas-expulsion delay times for all models are set to 
td = 0.6 Myr. The details of the models are given in Table 
1 . 


4 RESULTS 

In this section we present the results from our numerical sim¬ 
ulations. The mass-loss in the beginning of the evolution is 
dominated by gas expulsion and early mass-loss associated 
with stellar evolution of massive stars, such that clusters 
lose about 30% of their mass within the first 100 Myr by 


early stellar evolution. Moreover, in the framework of our 
model, the evolution of a star cluster is affected by several 
parameters related to the details of the gas expulsion pro¬ 
cess (like e.g., star formation efficiency, e , and gas expulsion 
time scale, tge ) as well as the strength of the external tidal 
field, quantified by the ratio of the initial half-mass radius of 
a star cluster to its tidal radius, r^/rt. In addition, the pres¬ 
ence of primordial mass segregation significantly affects the 
global dynamical evolution of star clusters. In tidally lim¬ 
ited clusters, primordial mass segregation leads to a stronger 
expansion, and hence a larger flow of mass over the tidal 
boundary llBaneriee fc KrouDall2014 1. It may therefore help 
to dissolve them more rapidly. Tidally-underfilling clusters, 
however, can survive this early expansion (Vesperini et al. 
2009). The degree of primordial mass segregation and the 
strength of the tidal field of the host galaxy are therefore a 
crucial parameter in the early evolution of GCs. We therefore 
compare the evolution of the MF of the clusters by changing 
the strength of the external field and degree of initial mass 
segregation. Table 1 gives an overview of the simulations 
performed. 


4.1 The influence of residual-gas expulsion on the 
MF of a primordially not mass-segregated 
cluster 

We first start with one model without primordial segrega¬ 
tion (’rhl0-S0-sfe25’). In order to see the maximum effect 
of residual-gas expulsion, we set the initial half-mass radius 
of the modeled cluster to Vh = 10 pc, and therefore it is a 
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Figure 1. Evolution of the global stellar mass function (for main-sequence stars) inside the tidal radius is computed at each snapshot 
for a tidally filling cluster without primordial mass segregation. The star formation efficiency is set to e = 25%. The dashed vertical 
lines show the positions of m = 0.5Mq and m = 0.85Mq. The values of the best-fitting slopes in the intervals, m ^ 0.5Mq (cti), 
0.5Mq ^ m ^ 0.85 Mq ( 02 ), and for m is O.SSMq ( 03 ) are written in the plots. As can be seen, the effect of gas expulsion on the mass 
function is negligible for initially not mass-segregated clusters. 


tidally-filling cluster (r^i/n = 0.18). We also set e = 25%, 
because the lower the SFE, the stronger is the effect of gas 
expulsion (GE). We plot the mass function at different snap¬ 
shots to see how it evolves with time in the intermediate- 
mass (0.50 ^ ^ 0.85) and high-mass (m 0.85Mq) 

part evolve with time (Fig.l). We also fitted the MF of 
main-sequence stars with masses in the three different mass 
ranges of m < O.hOM© (oi), 0.50 ^ ^ 0.85 (02), and 

m ^ 0.85Mq ( 0:3) with a power law function. 

As can be seen, the mass function of the initially 
non-segregated cluster does not change with time. Just a 
very slight decrease of the slope due to two-body relax¬ 
ation driven evaporation of stars occurs in the mass range 
0.50 Si ^ 0.85. This implies that the shape of the MF 
stays more or less the same, i.e. stars of all stellar masses are 
equally affected by residual gas-expulsion. In Figure [5] (see 
model ’rhl0-S0-sfe25’) we also show the time evolution of 
02- It can be seen that the MF flattening occurs on a much 
longer timescale for the non-mass-segregated system. 


4.2 The influence of residual-gas expulsion on the 
MF of primordially mass-segregated clusters 

We assigned a quite extreme degree of segregation of S = 0.9 
for all primordially segregated models, in order to show the 
maximum influence mass segregation can realistically have. 


4-2.1 Tidally-underfiUing clusters 

Because of stellar evolution post gas-expulsion clusters expe¬ 
rience a further expansion such that clusters with a higher 
degree of primordial mass segregation experience a larger 
jump in the half-mass radius within the first 100 Myr of 
evolution. Almost 30% of the initial mass is lost owing to 
stellar evolution on a timescale of 100 Myr. Fig. [2] depicts 
the effect of gas expulsion on the MF of a tidally-underfilling 
cluster with an initial half-mass radius of Vho = 1 pc and a 
filling factor of rh/rt = 0.02. As can be seen, the MF in 
the range m ^ 0.5Mq still resembles the IMF and does not 
show a strong change of the slope, while in the mass range 
m ^ O.5M0 there is a significant change in the slope. 

In fact, if the tidal field is weak (i.e. rh/rt = 0.02, see 
Table 1), even a model with a SFE as low as 25% keeps the 
majority of its mass and undergos very little mass loss due 
to the tidal interaction. According to the final mass demon¬ 
strated in Table 1, a very small fraction of mass is lost due to 
dynamical evolution caused by the gas expulsion process. In 
other words, the initially tidally-underfilling models experi¬ 
ence a weak tidal field, and hence keep most of their stars 
and retain their IMF slope in the mass range m 0.5Mq. 

Fig.|3]shows different snapshots of the mass function of 
a marginally tidally-underfilling cluster [rh/rt = 0.09 with 
r^o = 5 pc) with a primordial mass segregation with S=0.9 
and a star formation efficiency of e = 25%. As can be seen, 
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Figure 2. Evolution of the global stellar mass function within the tidal radius of an initially tidally underfilling cluster with primordial 
mass segregation (3=0.9). The star formation efficiency is set to e = 25%. The initial half-mass radius of the modeled cluster is r/jQ = 1 
pc. The dashed vertical lines shows the position of m = O. 5 OM 0 and m = O.85M0. The slope of the best fit line in this interval, 02 and 
for m ^ O.85M0, 03 , are written in the plots. As can be seen the slopes in these mass ranges do not change away from the initial form 
because the high-mass stars which are located in the inner part of such a compact cluster are less affected by the tidal interaction. 


gas expulsion significantly changes the slope of the mass 
function only in the low-mass part (m ^ O.5OM0). The 
mass function at the high-mass end does not change, while 
the slope in the intermediate mass range, O.5OM0 ^ m ^ 
O.85M0, changes from =2.3 to 2. Note that this model 
is not an extreme case of tidally-underfilling systems, but it 
is an intermediate case between the two regimes of tidally- 
filling and tidally-underfilling clusters. Figure [5] shows the 
time evolution of 02 for this model. 


2..8..S Tidally-filling cluster 

For clusters initially tidally-filling their Roche lobes (i.e. 
models with Vh/rt 0.1) and primordially mass segregated, 
a substantial fraction of the initially bound mass is lost sim¬ 
ply due to stellar evolution and expansion of the cluster 
driven by gas expulsion. Our simulations show that a given 
amount of mass loss leads to a stronger MF flattening than 
for tidally-underfiling clusters. The strong tidal fields lift 
even the clusters with a large SFE away from the initial 
IMF slope. 

Fig. |4] depicts the evolution of the slope of the mass 
function over the whole range of stellar mass at different 


snapshots for a tidally-filling, primordially mass-segregated 
cluster with rno = 10 pc, i.e. model ’rhl0-S0.9-sfe25’. The 
slope was determined from a fit to the distribution of stars 
in two different mass ranges, O.5OM0 ^ m ^ 0.85Mq (0:2), 
and m O.85M0 (03). Fig. 0] confirms that the slope of 
the mass function for tidally-filling clusters changes in the 
intermediate-mass range as well as in the low mass part, 
by getting shallower. Figure [S] shows the evolution of the 
MF-slope in the intermediate-mass range for all computed 
models. 

Comparing the last two models in Table 1 and Fig. O 
it can be seen that for initially segregated clusters, as the 
SFE decreases, the slope of the MF in the low-mass part 
decreases. This is easy to understand: The lower the SFE 
gets, the stronger is the dynamical effect of the gas expulsion 
process. 

The oscillatory trend of 02 in Fig. [5l particularly no¬ 
table for the mass-segregated tidally-filling clusters, are re¬ 
lated to the oscillatory trend of the mean velocity disper¬ 
sion (Fig. [9]) of the outer layers which are dominated by the 
low-mass population of stars (see Sec. 5 for more details). 
According to Fig. [6] after an early expansion driven by gas 
expulsion, the half-mass radius increases by a factor of 2 
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Figure 3. The evolution of the global stellar mass function within the time evolving tidal radius of an initially marginally tidally 
underfilling cluster with primordial mass segregation (S=0.9). The star formation efficiency is set to e = 25%. The initial half-mass 
radius of the modeled cluster is rf^Q = 5 pc (r;,/rt=0.09). The dashed vertical lines show the position of m = 0.5Mq and m = O.85M0. 
The slope of the best fit line in this interval, 02 (that changes from 02 = 2.27 at T=0 to 02 = 1.97 at T=100 Myr), and for m O.SMq, 
03 , are written in the plots. 


(i.e., Vh =20 pc), and the total mass within the tidal radius 
decreases to M = 20000Mq at T = 100 Myr, so the dynam¬ 
ical time scale of the cluster is equal to tcroaa — 30 Myr. 
This value is nearly equal to the period of the oscillations 
as can be seen in Figs. ISII^and llUl 


4.3 Comparison with observation: q — c plane 

Fig. [7] compares the influence of gas expulsion on the MFs 
of computed models, with the eye-ball fit to the observed 
data from De Marchi, Paresce & Pulone (2007, hereinafter 
MPP ) with the functional form of q;(c) = -1-2.5, supple¬ 

mented with additional global MF slopes taken from Paust 
et al. (2010). The canonical Kroupa IMF has a = 1.72 
over this mass-range (0.3 — O.SMq). The models with gas 
expulsion agree well with the general trend that less concen¬ 
trated elusters typically show a shallower MF. The non-mass 
segregated model (filled circle), and a test model without 
gas expulsion (filled square), are not able to reproduce the 
observed trend confirming the conclusion by Leigh et al. 
(2013). These results confir m, by direct N-body simulations, 
the co nclusions reached by iMarks. Kroupa. fc Baumgardd 
(l2008l) based on simplified models. In order to justify our 


conclusion we calculated five additional simulations without 
gas expulsion and primordial mass segregation for differ¬ 
ent values of the concentration parameters in the range of 
c = 0.85 — 2.65. The results are demonstrated as triangles 
in Fig. [71 As can be seen these models keep their initial 
MF-slope and independent of the value of the concentration 
parameters all lie above the observed data points. 


From Fig. [5l we see that the change in 012 over time is 
small/negligible for tidally under-filled clusters and is more 
noticeable for tidally marginally-filled and filled clusters. 
Note that the declining and oscillating trend of 02 for the 
tidally filling clusters in Fig. [S] would continue well beyond 
100 Myr and will moderately change their positions in the 
c — a diagram, they will still be consistent with the MPP 
trend since the latter has some scatter anyway. This is be¬ 
cause the estimated two-body relaxation times for all models 
with an altered MF are larger than one-third of the Hubble 
time at 100 Myr. Therefore, the 02 values will likely decline 
further due to unbound escaping stars over the intended 
mass range still remaining within rt- 
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Figure 4. The same as Fig. [3] but for a tidally-filling cluster with tjjq =10 pc. The slope of the MF in the mass range m = 0.50 Mq 
and m = O.SSMq significantly changes away from the canonical value of 2.30. The MF shows a nonlinear trend for m Sj O.SOMq, below 
which it becomes much shallower than the canonical IMF. An explanation for this specific form is represented in Sec. 5 


4.4 Implications for the very young MW 

Our simulations show that only in the case where a cluster 
ends up with a very large rhjrt, will the low mass stars leave 
the cluster. Note that all clusters in Table 1 which have a 
significantly altered MF due to gas expulsion end up with 
final 3D half-mass radii of about ~21 pc. This is compatible 
with the typical size of remote halo globular clusters such as 
Palomar 4 and Palomar 14, but they are still larger than the 
size of the clusters in the MPP plot with a median half-mass 
radius around < >= 8.9 ± 5.9 pc. 

Although the present-day radii of the MPP clusters do 
not fit a tidally-filling scenario, this does not mean that the 
particular clusters that have a damaged MF were not tidally 
filling during the fist 10-100 Myr of their life when the expul¬ 
sion of residual gas caused violent relaxation and expansion. 
The present-day population II halo GCs which formed with 
a higher metallicity, and thus typically after the GCs with a 
lower metallicity, most probably formed in an environment 
which had larger fluctuations in the tidal field (Marks & 
Kroupa 2010, their fig. 12). 

Also, during the first few Gyr the MW did not exist 
in the way we know it today. Observations have revealed 
that the size and the mass content of galaxies change sig¬ 
nificantly with redshift such that the sizes of the galaxies 
at high redshifts are smaller in comparison with galaxies of 


similar mass in the local universe 

e.e.. iFranx et al. 

2008; 

Williams et al. 2010l: Mosleh et al.l 

201ll: iLaw et al. 

2012; 

Mosleh. Williams & Franx||2013|). I.e., the fluctuation in the 


pre-MW potential grew with time as it began assembling 
during the first few Gyr and it is only during the past 8-10 
Gyr that the MW finally settled down in the form we know 
it today. So, when the GCs were relatively young, their envi¬ 
ronment would have been much more dramatic. Some, but 
not all, of these GCs would have formed in a strong lo¬ 
cal tidal field such that they lost a significant fraction of 
their stellar population as a result of the expulsion of their 
residual gas. For example, a very young GC with a mass 
of Mgc = 10^ Mq would have formed with an embedded 
half-mass radius of about 0.8 pc (assuming the radius-mass 
relation of Marks & Kroupa 2012). Expulsion of 67 per cent 
of unused gas would lead to an expansion to a half-mass 
radius of about 3 pc, and further expansion through stellar 
evolution would continue. If the IMF was top-heavy (Marks 
et al. 2012), then the expansion would be larger still. If this 
young GC was mass-segregated and immersed in a tidal field 
from a nearby extremely massive molecular cloud with a 
mass of, for example, Mmc = 10® Mq , which is in the pro¬ 
cess of forming another GC, then the tidal radius would be 
rt ~ D (Moc/MMcf^^ = 23 pc such that this GC would 
have had Vh/rt > 0.13 if it was located D = 50 pc away 




























































10 Haghi et al. 


3 


2.5 


V 

E 

V 
© 

5 

u, 1'5 

o 

!?’ 

1 


0.5 

0 20 40 60 80 100 

Time [ Myr ] 

Figure 5. The time evolution of the MF-slope, 02 , for main se¬ 
quence stars within the tidal radius in the mass range 0.50Mq 
M ^ O.85M0. The initial number of stars in all models is 
N=170000 (corresponding to M = 10®Mq) with an input Kroupa 
canonical IMF with the slope of 02 = 2.3. The initial half-mass 
radius, star formation efficiency, and degree of primordial segre¬ 
gation (i.e., S = 0 and S = 0.9) are indicated in the model name 
in the legend. As can be seen only tidally filling clusters with 
primordial mass segregation show a significant flattening in the 
slope of the MF. In the primordial mass segregated models, the 
lower SFEs lead to a shallower MF-slope. 
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from the cloud. Such events and configurations would have 
persisted during the first 0.5-1 Gyr of the evolution of the 
proto-MW and the above case would exist for 10-100 Myr 
given that the just-hatched GC would be approximately co¬ 
moving with the evolving nearby molecular clouds during 
half an orbit through the forming proto-MW. Present-day 
instances of such processes and events are evident in the An¬ 
tennae galaxies where massive star cluster complexes con¬ 
taining hundreds of massive cl usters within region with a 
diameter o f 300 pc are formin g dKroupal 1998 Wilson et al.l 
I 2 OO 0 I . I 2 OO 3 I : iBriins et al.ll201l] : I Johnson et al.l 20151) . 

Since the young GCs did not form in a smooth environ¬ 
ment, they must have been exposed to strongly fluctuating 
tidal fields in the very early pre-MW potential. ”R136” is a 
good example: it is the most massive compact cluster in a 
whole cluster complex. There are many lesser clusters that 
formed in the ”R136” star-burst region. The ’’Arches” clus¬ 
ter is in fact in a very strong tidal field, and most probably 
also did not form alone. However, these cases are almost two 
or three orders of magnitude lesser in mass scale than what 
was occurring in the very young pre-MW. So, the damaged 
GCs (i.e., the ones with low c and depleted MF) may well 
have been the lesser clusters in a star-burst region orders of 
magnitude more intense than anything in the present-day 
local Universe. 

It should be noted that the correlation is not a pure 
time-sequence. But the more damaged clusters are typically 
more metal rich (Marks & Kroupa 2010). This means they 
typically formed a few hundred Myr later than the first, 
metal poorest clusters which are also least damaged and 
would have formed in an overall smoother potential (weaker 
tidal fields) and more compact due to the lower metallicity. 

The higher—[F’e/IT] clusters are further inside the MW, 


rf,Q=10 pc, e=25%, S=90% 



Time [Myr] 

rf,g=10 pc, 0=25%, S=90% 



Figure 6. The evolution of total mass (upper panel) and half¬ 
mass radius (lower panel) with time for simulated star cluster 
”rhl0-S0.9-sfe25”. The cluster loses about 80% of its mass by 
early rapid gas expulsion and stellar evolution. 


which means they formed slightly later and under more vi¬ 
olently varying tidal fields, leading to the MPP correlation. 
Note that this correlation only means that at low c (typically 
larger [Fe/H]) there is a wider spread of alpha values than at 
high c (typically lower [Fe/H]). The higher-metallicity clus¬ 
ters typically but not exclusively formed more tidally filling 
due to larger mass concentrations nearby them, and they 
formed somewhat deeper in the proto-MW potential. 

Other mechanisms may contribute to the MPP trend. 
For example, by assuming a high retention fraction of 
stellar-mass BHs after supernovae which prevent the clus- 
ters from going into core collapse for nearly a Hubble time, 
iLiitzgendorf, Baumgardt, fc KruiissenI (l2013l) explained the 
MPP relation without invoking primordial mass segregation 
and gas expulsion. They found that a cluster hosting a high 
number of stellar-mass black holes increases the escape rate 
of high-mass stars from a star cluster, implying that the rel¬ 
ative depletion of the mass function at the low-mass end 
proceeds less rapidly. 

Two-body relaxation plus a tidal field, and any pro¬ 
cess that can modulate the cluster expansion according to 
concentration (more concentrated clusters expand less) can 
produce a MPP like trend. Furthermore, the observed metal¬ 
licity trend of the MPP clusters (more damaged GCs have 
typically higher [Fe/H]) needs to be understood in terms 
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Figure 7. The slope o of the MF at the end of the N-body in¬ 
tegrations (T = 100 Myr) in the mass-range 0.3 - 0.8 Mq as a 
function of the cluster concentration, c = logjQ(rt/rc), for all cal¬ 
culated models listed in Table 1 (circles with central dots). The 
solid line is an eye-ball fit to the observ ed data (cros s es) fro m De 
Marchi, Paresce &; Pulone (2007) and IPaust et al.l ll201Cll . The 
filled square marks the primordially segregated modeled cluster 
without gas eiprdsion,”rhl0-S0.9-NoGE”, the filled circle repre¬ 
sents the non-segregated model ”rhl0-S0-sfe25”, and the filled tri¬ 
angles represent simulations without gas expulsion and primordial 
mass segregation. 


of the stellar-mass BH scenario. In the long term, neutron 
stars (NSs) will also influence the cluster’s expansion. It is to 
be noted that the BHs and NSs influence mostly the central 
part of the cluster, the expansion (and contraction later) 
becomes less prominent as one goes outwards. In any case, 
BHs (and NSs) can keep the cluster from core collapsing for 
a long time. Their effect in terms of metallicity can be more 
subtle. In fact, two opposite trends seem to be there in the 
a — [Fe/H] correlation: 

• With lower metalicity, BHs are more massive and would 
expand the cluster more initially and that tends to produce a 
trend which is opposite to the MPP’s a—\Fe/H] correlation 
(Marks & Kroupa 2010). 

• On the other hand, more massive BHs eject each other 
more efhciently and, in the long run, will be more depleted 
which is favorable for the metallicity trend. 

The asymptotic behaviour of stellar remnants in star clus¬ 
ters and their effects are still largely unexplored. However, 
straightforward canonical stellar dynamical models have 
been implying that the MPP relation is unexpected (Leigh et 
al. 2013) unless gas expulsion and mass-segregation at birth 
with a strongly varying tidal field is taken into account. 


5 EVIDENCE FOR PRIMORDIAL MASS 
SEGREGATION 

As we have shown in a strong-enough tidal environment, 
a strongly flattened slope of the stellar MF in the low- 
(m ^ O.SMq) and intermediate-mass (~ 0.50 — O.SSM©) 
regime results after the rapid expansion driven by residual- 
gas removal and massive star evolution. As can be seen in 
Figs. l2l4l such a flattening is the exclusive feature that occurs 


only for initially mass-segregated clusters for both regimes, 
tidally filling and tidally underfilling. Therefore the imprint 
of residual gas expulsion, as direct evidence of primordial 
mass segregation, might be visible in the present-day MF. 
Therefore, this can be interpreted as an observational con¬ 
sequence of gas expulsion o n the MF in the low mass end 
confir ming the conclusions of iMarks. Kroupa. fc Baumgardtl 

ll2008l) . 

It should be noted that, apart from the particular shape 
(the dip near ~ 0.2Mq) of the MF in the mass range m ^ 
0.5Mq, a cluster which follows the MF shown in, like e.g.. 
Fig. 4 {{rh/rt)o = 0.18, with S = 0.9) would end up with 
a mass-to-light ratio smaller by a factor of 0.6 (due to the 
depletion of low-mass stars and remnants) compared to the 
case of Fig 1 {{rh/rt)o = 0.18, with S = 0.0). Of course, 
this conclusion depends on the way we set up the initially 
segregated clusters. If we initialize clusters such that there is 
not such a large degree of mass segregation among the low- 
mass stars, but only the high mass stars are near the cluster 
centre, the results will be changed. In this section we detail 
our calculations to find an explanation for the particular 
shape of the mass function in the mass range m ^ 0.5M©. 

Fig. [8] (right panel) shows the evolution of the 3D half¬ 
mass radius of 6 different mass-bins within the range 0.08 
to 0.8 Mq for a computed model with primordial mass seg¬ 
regation. For comparison purpose we compute an identi¬ 
cal model without primordial mass-segregation (left panel). 
Both models are tidally filling with an initial concentration 
of {rh/rt)o = 0.18. In the not mass-segregated model, all 
stars of different mass-bins are equally affected by early 
gas removal, having the same half-mass radius and they re- 
virialize in 20 Myr. But, in the case of the primordially mass 
segregated model, the outer parts which are filled by stars in 
the mass range m ^ 0.2Mq fall back into the inner regions, 
while the next inner shells containing the stars with higher 
masses (m 0.2Mq ) move outwards enhancing the escape 
rate of stars in this range from the cluster. 

Fig. 8 shows why a minimum can be seen in the stellar 
mass function near 0.2Mq with a rise towards lower masses 
(Figs OH): In a fully mass segregated and virialised initial 
cluster the least massive stars (?( O.IM©) occupy the least 
bound orbits at largest radii such that they can only fall to¬ 
wards the cluster centre. They do this despite the expulsion 
of residual gas. More massive stars are lost from the clus¬ 
ter due to gas expulsion, while the least massive stars are 
preferentially accumulated. The initial separation of half¬ 
mass radii of the different mass-bins vanishes as a result of 
phase mixing and two-body relaxation and all layers tend 
to have the same steady half-mass radius long-term after 
re-virilization. 

Fig. [9] shows the evolution of the mean velocity disper¬ 
sion within the tidal radius of the 6 mass-bins over 0.08 - 
0.80 M© for the same models as shown in Fig.jS] In the non- 
segregated model, the mean velocity dispersion of all stars 
with different masses equally decreases because of the rapid 
expansion of the star cluster caused by rapid gas removal. 
In the case of the primordially mass segregated model, and 
as explained above, the least massive stars begin to fall to¬ 
wards the cluster centre. Their velocity dispersion increases 
while the more massive stars are partially lost from the clus¬ 
ters. The least massive stars retain their low-energy orbits 
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after re-virialization thereby oscillating through the cluster. 
Their velocity dispersion therewith oscillates. 

The corresponding evolution of mean radial velocity 0 
within the tidal radius is demonstrated in Fig. llOl As can be 
seen all the shells in the not mass-segregated cluster expand 
equally and after 20 Myr they reach a state of equilibrium. 
The initial increase in the radial velocity is due to dilution 
of the gravitational potential well driven by the early gas 
expulsion. In fact, in a virialized system, the stars are bound 
in different orbits and the mean radial velocity of all stars is 
zero. But, the early ejection of gaseous mass causes a cluster 
to expand and leads to a positive mean radial velocity of 
the stars. For initially mass-segregated clusters, the mass 
lost due to the evolution of the massive stars is removed 
preferentially from the cluster inner region, and the early 
expansion of the cluster is stronger. 

In order to quantify the effect of gas expulsion, we com¬ 
pare the evolution of the mean radial velocity in different 
mass bins for a model without gas expulsion (top panel of 
Fig, mil with the gas expulsion model (bottom panel of Fig. 
El) which is otherwise equal to the case shown in the top 
panel of Fig. El Since the least massive stars (that are re¬ 
siding in the outer parts of the cluster) have the smallest 
velocity dispersion, they are least affected by gas expulsion. 
The lower mass bins gain a larger mean radial velocity due 
to gas expulsion leading to an expansion. So, they are more 
prone to escape from the cluster. 

The evolution of the half-mass radii of different mass- 
bins for both models as a function of time is plotted in Fig. 
El Note that in the gas expulsion model, the initial sepa¬ 
ration of half-mass radii of different mass-bins is reduced as 
a result of violent relaxation driven by rapid residual gas- 
expulsion. But, the model without gas expulsion preserves 
this initial separation and is still mass segregated at T = 100 
Myr. 


6 CONCLUSION 

We have performed a series of N-body computations study¬ 
ing the impact of expulsion of residual gas on the stellar 
masses function slope in star clusters for a wide range of 
stellar mass, varying the strength of the external tidal field 
in order to model the large spatial differences of the tidal 
held expected in the varying proto-MW (Marks & Kroupa 
2010 ). 

Gas expulsion leads to the loss of a substantial frac¬ 
tion of the initially bound mass. A tidal held and mass loss 
by stellar evolution helps increasing this unbound fraction. 
Comparing the hnal masses in the hrst and last two mod¬ 
els in Table 1, one can see that the mass depletion due to 
gas expulsion can be comparable or even larger than that 
due to stellar evolution, particularly for clusters subject to 
a stronger tidal held modelled here as tidally hlling clusters. 

We have shown that the MFs of non-segregated clusters 
still resemble their IMFs and that they do not show a strong 
change of the slope over the whole range of stellar mass over 
100 Myr of evolution, because stars of all stellar masses are 
equally affected by gas expulsion. 

® The mean radial velocity is defined as < v.f > 


Inatially mass segregated without gas expulsion 



Time [ Myr ] 

Initially mass segregated with gas expulsion 



Figure 11. The evolution of the mean radial velocity within 
the tidal radius of stars in different mass-bins for two computed 
initially mass segregated tidally-filling models, ’rhl0-S0.9-NoGE’ 
(top: without gas expulsion) and ’rhl0-S0-sfe25’ (bottom: with 
gas expulsion). The models are the same, except that one is w/o 
gas expulsion. The least massive stars, that initially reside in the 
outer parts of the cluster, are least affected by gas expulsion since 
they have small velocities. 


In the case of primordially mass segregated clusters, the 
models divide into two different regimes, depending on the 
strength of the external tidal field: 

• Tidally-underfilling models with rh/rt values smaller 
than 0.1 (i.e. when the cluster size is small at first) which 
keep their MF in the mass range m ^ O.SOM© unchanged. 

• Tidally-filling models with rh/rt values larger than 0.1, 
where a large fraction of their low-mass stars are lost and 
hence they evolve to a comparatively shallow MF. 

Clusters with ru/rt ^ 0.1 approximately preserve 
their MF at the high-mass end (i.e. m ^ 0.85Mq) apart 
from changes due to stellar evolution, and receive only a 
very small change in the intermediate stellar mass range 
(0.50Mq ^ m ^ O.85M0) even if the SFE is low (SFE 
^ 0.25%). But the MF in the low-mass part (m ^ O.SMq) 
changes significantly, because primordial mass segregation 
leads to a larger impact of stellar evolution on cluster ex¬ 
pansion and enhances the preferential loss of less massive 
stars due to two-body relaxation after post-gas expulsion 
re-virialisation. 
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Figure 8 . The evolution of the half-mass radii of 6 mass-bins chosen within the range 0.08 to 0.80 -ATq, for two computed tidally-filling 
models, ’rhl0-S0-sfe25’ (left: without primordial mass segregation) and ’rhl0-S0.9-sfe25’ (right: with primordial mass segregation). The 
half-mass radius of the cluster is also plotted for comparison. Note that in the non-segregated model, stars with all stellar masses are 
distributed with the same half-mass radius being equally affected by gas expulsion. In the case of the primordially segregated model, the 
stars in the mass range m ^ O.ISMq move inward during the first 10 Myr and are hence less prone to escape from the cluster, while the 
next inner shells containing the stars with higher mass (O.ISMq ^ m ^ 0.25Mq ) move outwards after an initial brief contraction phase 
strengthening the possibility of escape from the cluster. This is reflected in the stellar mass function of the cluster (Figs. 3-5) as forming 
a minimum at m = 0.2Mq. Note that the initial separation of half-mass radii of the different mass-bins is reduced as a result of violent 
relaxation driven by rapid residual gas-expulsion and all layers tend with increasing time towards the same steady half-mass radius after 
re-virilization. 
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Figure 9. Left: The evolution of the velocity dispersion of stars within the tidal radius for different mass-bins as in Fig. [8] Left: For a 
non-segregated model the global velocity dispersion of all stars in all mass-bins equally decreases with time (no difference for different 
mass-bins) owing to a rapid expansion with timescale tqe driven by gas expulsion. After a few crossing times the cluster reaches its 
steady equilibrium velocity dispersion. Right: the same as the left panel but for the primordially mass segregated model. Note that 
the inner shells (i.e., the higher mass bins) initially have a larger velocity dispersion consistent with the virial equilibrium state. After 
gas-expulsion and corresponding dilution of the gravitational potential well, owing to the more rapid expansion of the inner shells (with 
an initially higher velocity dispersion), the order of the layers is inverted in terms of their velocity dispersion. The oscillation arises from 
the outer low-mass stars orbiting through the expanded cluster which has a two-body relaxation time of 10 Gyr (Table 1). 


On the other hand, primordially mass-segregated 
tidally-filling clusters with Vh/rt ^ 0.1 show a strongly de¬ 
pleted mass function in the intermediate stellar mass range 
(O.5OM0 ^ m ^ O.85M0) as well as in the low-mass part. 
Also, the slope of the MF at the high-mass end changes 
slightly. 

The mass functions in both cases are quite differ¬ 
ent from the non-segregated models, especially in the 
intermediate-mass and low-mass part: moving towards lower 


stellar mass, the number of stars per mass interval at hrst de¬ 
creases and then starts to increase around m ~ 0.2Mq (i.e. 
it is not possible to fit the MF with one power-law function). 
If such a behavior in the MF-slope of the low-mass part is 
observed, this would constitute direet evidence for mass seg¬ 
regation, gas expulsion and a strong tidal field at birth. 

Concerning the dip at the low-mass end of the stellar 
mass function of fully mass segregated initial clusters, our 
explanation is that the least massive stars are essentially on 
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Non-mass segregated Mass segregated 



belongs to the tidally-filling non-mass segregated cluster with rho =10 pc and the right one shows the results of a similar cluster but 
with primordial mass segregation. 


radial orbits, since they cannot be moving outwards much 
given that the cluster is in virial equilibrium. The least mas¬ 
sive stars must have the smallest velocity dispersion. So, 
they probably fall in, or at least, they are least affected by 
gas expulsion since they have small velocities, while the more 
massive stars are lost from the cluster due to gas expulsion. 

It should be noted that if the outermost stars predomi¬ 
nantly move inwards, that is more of an artefact of how the 
mass-segregated cluster is set up. However, even without the 
infalling nature, a similar feature would have appeared since 
the outermost (and lowest mass for a fully mass-segregated 
cluster) stars are the slowest. 

If gas expulsion is a common process in the early stages 
of star clusters, then by direct measurement of the slope 
of the mass function at low stellar mass of very young 
massive clusters 50 Myr), one can constrain the start¬ 
ing condition of the clusters such as primordial mass seg- 
regation. These results confirm the c onclusions reached by 
iMarks. Kroupa. fc Baumgardtl ll2008l b Moreover, the shape 
of the MF can constrain the birth place or the birth size of 
the cluster, because the tidally-underfilling clusters have a 
falling and then a rising trend in the MF-slope in the low- 
mass end if they are born initially segregated. 

The models with gas expulsion agree with the general 
trend that less concentrated clusters show a flattened MF. 
We note that the MF-slopes and c values are determined in 
the models presented here after 100 Myr of evolution, but, 
since the two-body relaxation time for all models with an 
altered MF are larger than one-third of the Hubble time 
(Table 1), the position of the clusters in the c — a plane will 
not significantly change over 10 Gyr. 

Our models constitute a possible solution to the flat¬ 
tened mass function of the two remote halo GCs, Pal 4 and 
14 (Zonoozi et ah, in preparation). In our previous publica¬ 
tions (Zonoozi et al. 2011, 2014) we have shown that, while 
two-body relaxation does have an effect on the cluster’s stel¬ 
lar mass function, the effect from two-body relaxation is too 
weak such that the shallow present-day mass function slope 
measured by Jordi et al. (2009) and Frank et al. (2012) can¬ 
not be reproduced with models starting from a canonical 
Kroupa IMF (2001). We have demonstrated that the ob¬ 


served parameters of globular clusters Pal 4 and Pal 14 (i.e. 
their observed mass functions, mass segregation, half light 
radii and masses) can, in principle, be reproduced either by 
assuming a very high primordial mass segregation, or by 
assuming an IMF already depleted in low-mass stars, com¬ 
pared to the canonical Kroupa IMF. We showed that such an 
initially flattened mass function can be explained by the loss 
of a good fraction of low-mass stars from an initially mass- 
segregated embedded cluster across the tidal radius during a 
violent early residual-gas expulsion phase in agreement with 
the general findings by Marks & Kroupa (2010). 

Therefore, the flattened IMF of Pal 4 and 14 we have un¬ 
covered strengthens the idea that these clusters formed with 
a high degree of primordial mass segregation in a stronger 
tidal field rather than in isolation. I.e. if Pal 4 and 14 have 
always been on a circular orbit with the present day Galac- 
tocentric distance, it is hard to imagine that the tidal field 
would have played a role during the response of the cluster 
to gas expulsion. This may have been the case if either the 
Gala ctic tidal field has since then evolved dMarks &: Kroup^ 
l2010lf . or if the cluster was born within the inner part of our 
Galaxy, and then moved to a larger Galactocentric distance 
on an eccentric orbit, which is more likely than a circular or¬ 
bit. Or, the cluster birth site was a now detached/disrupted 
tidal dwarf galaxy (Pawlowski et al 2011). Evidence for this 
scenario would strengthen if these clusters are found to or¬ 
bit within the vast polar structure (VPOS, Pawlowski & 
Kroupa 2014). 


Thus the two most important new results we arrive at 
are (i) that the observed correlation between concentration 
and the MF-slope can be accounted for excellently (Fig. 7) 
and (ii) that the MF of stars ought to have a minimum at 
about O.2M0 (at an age up to ~ 100 Myr), if GCs formed 
mass-segregated and suffered gas expulsion in a strongly 
varying proto-Galactic pote ntial. This confirms the previ¬ 
ous co nclusi ons arrived at by Marks. Kroupa. fc Baumgardtl 
(l2008ll and iMarks fc Kroupal ll2010ll . respectively, yielding 
an important if not essential window towards studying the 
very early phase of MW assembly. The approach taken here 
(gas expulsion and tidal field) yields otherwise unobtainable 
constraints on the very early conditions when the proto-MW 
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Inatially mass segregated without gas expulsion 



Time [ Myr ] 
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Figure 12. The evolution of the half-mass radii of 6 mass-bins 
chosen within the range 0.08 to 0.80 Mq, for two computed 
tidally-filling models as in Fig. 1111 The half-mass radius of the 
cluster is also plotted for comparison. 


was forming. That is, with this approach, we can obtain con¬ 
straints on a time-resolution of less than a few hundred Myr 
and on the fluctuations of the very early MW potential. 

We emphasize that massive embedded clusters with 
masses larger than about 10"* Mq revirialise rapidly and 
within a few Myr, they expand by a factor of 3-4 and loose 
a small fraction (at most 20 % in a solar neighborhood tidal 
field) of their stellar mass as a result of the blow-out within 
a thermal time-scale (with a gas velocity of 10 km/s) of 67 % 
of their original embedding gas mass (Brinkmann, Banerjee 
& Kroupa, in prep.; Banerjee & Kroupa 2014, 2015). Less 
massive clusters expand more, take a much longer time to re¬ 
virialise and are much more damaged losing more than 50 % 
of their stars by gas blow-out as shown by Kroupa, Aarseth 
& Hurley (2003). Thus, for those globular clusters, such as 
studied in this paper, which have a deficit of low-mass stars 
to be sufficiently damaged by gas blow-out, they must have 
been also subject to a strong local tidal field for the outer ex¬ 
panded regions of the cluster to be stripped. Only a fraction 
of all globular clusters suffered this fate (Marks & Kroupa 
2010). 

It should be pointed out that our conclusions in the 
present paper are based on two major assumptions, that 
is, primordial mass segregation and residual gas expulsion 
that are questioned by some authors dKruiissen et al.ll2012l : 


P arker. Dale, fc Ercolano 

[20121. I2OI5I: ISmith et ^ 


[20151 [Dale, Ercolano. fc Bonn^ 


201 if) . Concerning the debate 


whether gas expulsion is relevant in the formation of 
star clusters (Sec. 2), we note the high-precession space- 
astrometry of the GAIA mission is ideally suited to set¬ 
tle this argument, by allowing the kinematical fields around 
young clusters to be mapped accurately (Kroupa 2005). Sim¬ 
ulations with gas expulsion of clusters in the present-day 
MW need to be done to make predictions of the velocity field 
for GAIA, because each young cluster should be surrounded 
by a radially expanding population of stars of similar age as 
the cluster population. As predicted by Kroupa (2005) the 
expanding population forms a new moving group stemming 
from explosive residual gas expulsion that is related to the 
pre-gas expulsion velocity dispersion of stars in the embed¬ 
ded cluster, and thus to the mass of the re-virilised cluster. 
This new moving group (MGI according to Kroupa 2005) is 
present in addition to the classical moving group (MGII ac¬ 
cording to Kroupa 2005) stemming from secular evolution. 
Glustered star formation may thus also affect the kinematics 
and structural properties of galaxies. For example Kroupa 
(2002) has shown that the thick disk may be outcome of 
violent star formation within a thin disk. 
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